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1 Introduction. General equations 

The long wavelength excitations of a quark-gluon plasma are collective excitations which 
are described by nonlinear equations generalizing the classical Yang-Mills equations in 
the vacuum. Most studies have been limited so far to the weak field limit, where the 
equations become linear and the excitations reduce to abelian-like plasma wavesfll], very 
much similar to the electromagnetic waves in ordinary plasmas|0]. The purpose of this 
paper is to present new, truly non-abelian, plane-wave solutions that we have obtained 
recently. At leading order in the gauge coupling g, (we assume (? ^ 1 in the high 
temperature, deconfined, QCD plasma), the collective dynamics is entirely described by 
a set of equations for the gauge mean fields A^{x) which describe the long wavelength 
(A ~ i/gT) and low frequency {u ~ gT) excitations (T denotes the temperature) p, ^. 
(Throughout this work, the greek indices refer to Minkovski space, while the latin ones 
are color indices for the adjoint representation of the gauge group SU{N).) The equations 
satisfied by A^{x) are 

[D'',F^,ix)l=j!i{x), (1.1) 

where D^' = d>' + igAi^ix), At" = A^^T", and F>"' = [D'', D'']/{ig) = F^'^T''. In the right 
hand side, jj^ is the induced current which describes the response of the plasma to the 
color gauge fields A'^. (We do not consider here any external color source.) Its expression 
isi,i 

faix) =3u;lJ^v^W'Ax;v), (1.2) 

where Up = {2N + Ni)g'^T'^/18 is the plasma frequency, v^ = (1, v), v = k/k, k = |k|, 
and the integral / dQ runs over all the directions of the unit vector v. The functions 
W^ = W^T°- are defined as the solutions to 

[vD^^Wix^v]] = F^P{x)vp. (1.3) 

For retarded conditions (i.e., ^^(x) ^ as Xq ^ — oo) 



oo 



W^{x]v) = duUab{x,x-vu)F^^%x-vu)vp, (1.4) 



where U{x, y) is the parallel transporter along the straight line 7 joining x and y: 

U{x,y) = Pexp{-tg [ dz^A^iz)}. (1.5) 

The energy-momentum tensor of an arbitrary gauge field configuration in the plasma has 
been recently computedf^, with the following results for the energy density 

T'\x) = 1(E'^(x)-E-(x) + B«(x)-B'^(x)) + lu;lJ^W'A^;v)W'A^;v), (1.6) 



•iO 



and for the energy fiux density, or Poynting vector, S^ = T' 

S{x) = E«(x) X B«(x) + 1^1 J^^ W'^{x;v)W'^{x;v). (1.7) 

The non-abelian field strengths are, as usually, E^ = F*° and -B^ = —{l/2)e^^^Fl^. The 
first terms in the r.h.s. of eqs. (|1.6| )-( pr7| ) represent the standard Yang- Mills contributions; 
the terms depending on W2{x\ v) are related to the color polarizability of the plasma. In 
the absence of external sources, the energy conservation requires that 

9or°°(x) + diS\x) = 0. (1.8) 

It can be easily verified by using the equations of motion ( p,.l|) -( pr3D that this equation is 
indeed satisfied. 

2 Plane-wave solutions 

In this letter, we study particular plane-wave solutions which depend on x^ only through 
the variable z = p'^x^j,, where p^ = {00, p) is a fixed, time-like, four- vector {p^Pfj, = uj'^—p^ > 
0, p = IpI). Thus, we search for solutions of the form 

A^(x) =A^,{p-x), (2.1) 

in a gauge to be specified later. For the classical Yang-Mills equations in the vacuum, 
plane- wave solutions of this form have been investigated in Refs. 0, |^ . Another class of 
non-abelian plane- waves in the vacuum was considered by ColemanH]. We note that, in 



contrast to the vacuum case, the solutions of eqs. ( |1.1| ) for the high temperature plasma 
have direct physical relevance: they correspond to the collective color excitations of the 
QCD plasma. 

In deriving eqs. (^TT])-(|Op, we have assumed^ the gauge fields to be weak {A ^ T) 
and slowly varying {dA ~ gTA). For consistency, we therefore require that cj ~ p ~ gT 
and l^^l ^ T for the plane waves (|2.1| ). An important consequence of the Ansatz ( p.l|) is 
that the corresponding functions W^{x; v) are local and linear in the gauge fields: 

W^^ix; v) = - ^^(^) + ^vAa{z), (2.2) 

V ■ p 



as can be verified on eq. ( |1.3| ). Then, the induced current ( |1.2| ) may be easily computed, 
with the following results 



Pa{x) = - 3col aico/p) U^ - ^ ^°) , (2.3) 

and 

Ja(^) = -3.;J^{a(c./p)(^^ - 5^°)P - /3(^/p) ^r}, (2.4) 

where A^ = p • A-a (p = p/p), and Aj = A-a — P^a- The functions a{u) and /5(m) 



[u = uj/p) appear when evaluating the angular integral in (|1.2|) . We have {uj > p) 



Air 00 — V ■ p p^ p'^ 



(2.5) 



with 



a{u)^u(Q{u) -l), /5(n) = |(l - !^g(n)), (2.6) 



u^ 



and Q{u) is defined by 



QM-^ln^. (2.7) 

2 u— I 



For M > 1, the functions a{u) and I3{u) are positive. Remark that the current ( |2.3| )-( [27^ ) 
is not only covariantly conserved, [D^, j'^] = 0, as it should for the consistency of eq. ( |1.1| ), 
but it also verifies d^j'^ = p^Ja = 0. 



The expressions ( |2.3| )-( |27^ ) for the induced current are formally identical to those 
corresponding to an abelian plasma, that is, they involve only the gauge field polarization 
tensor, and no higher order irreducible amplitudes. This simplification represents an 
essential feature of the plane- wave fields that we are considering here. 

For later use, we evaluate here the polarization pieces of the energy density ( |1.6|) and 
of the Poynting vector ( |1.7|) for the plane- waves (|2.1| ) . We need the integrals {uj > p): 



Att {uj — v ■ p) 



"^'■p'p> + ^-jf^{5'^ -p'p^), 



pz. 



p^ 



and 



where 



and 
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yl yj yk 



An (tu — V ■ p) 



^ ' -p'p^p" -\ ' {d^^p' + d'-'^p' + d'^ p'') 



a(u) = 1 + 



p. 



u 



p. 



m2-1 



2Qiu) 



h{u) = Qiu) 



{2.1 



(2.9) 



(2.10) 



d{u) 



^u{Q{u) - l) - ^Q{u) 



c{u) + 3d{u) 



u 



u^ 



2 ^"' ' ' 2u 
1 - 3{u^-1){Q{u) - 1 



(2.11) 



with Q{u) given by ( p.7| ). For u > 1, the functions a{u) and b{u) are both positive. By 



using these results, together with eq. (2.2), we get 



T^'(x) 



1 



E"-E° 



UJ 
p2 



B"-B1 + i^P S^(^/^)M' - -<y + K-^lv)Al-Al 



UJ 



(2.12) 



and 



S(x) = E-^ X B'^ + ^u;2 ^|(^c(^/p) + 3^(^ 






P 



+ d{ujlp)J^,.Alv^2d{ujlp){Al 



^^°U^ 



ijj 



a •'a 



(2.13) 



In the rest of this letter, we work in the covariant gauge p^A^^ = 0. With the Ansatz 
A]), the gauge field equations ( |1.1| ) become then 



p^p^A^ + igp^ Ay, A 



A 



w-^ , ./~\. 



(2.14) 



(Throughout, the overdots indicate derivatives with respect to the argument of the func- 
tion, here z.) When contracted with p^, eq. ( p. 141) reduces to 



(■/t-j/ 5 w'H. U 



(2.15) 



(recall that p^Pu > and p^jf^ = 0). Therefore, we need only consider the simplified 
equation 



p^p'A'' - g' 



•y~\.y 5 •y~\. 5 •y~\. 



3' 



(2.16) 



together with the constraint (|2.15|) . In the present gauge, ujA^ = pA^, and the induced 
current (12. 31) -(12. 41") can be rewritten as 



Pa 
Ja 






O^ AT 



(2.17) 



where 



and 



»l^H^ 


P' 


Q - - 1 
V \p) 


02 - ^, ,2^^ f^ 


— 


u^ \p 



(2.18) 



(2.19) 



Since uj > phj assumption, the r.h.s.'s of eqs. ( |2.18| ) and ( p.l9| ) are both positive, so that 
Ql and Qt are real functions of uj/p. 

Any solution of the mean field equations (|2.15| )- (p.l9|) depends parametrically on the 
wave- vector p^, i.e. A'^ = A'^lp-x'jPu). It can be easily verified that A'^ is homogeneous in 
p'^ of degree zero, i.e., A^{\p ■ x; \p^) = A^{p ■ x;p^) for arbitrary constant A. Therefore, 
the solution depends only on three parameters, for instance, p'^/oo, with i = 1, 2, 3. 

At this stage, it is convenient to introduce three polarization vectors e'^(p;s), with 
s = 1, 2, 3 and p^e^(p; s) = 0. We choose the vectors e^(p; s = 1, 2) transverse to p, i.e.. 



e''{p; s) = (0, e,(p)), p ■ e^ip) = 0, e,(p) • e'^(p) = 5s 



(2.20) 



while 

e^(p;3) = ^=L=(p,^p). (2.21) 

The normahzation is such that e{p; s) ■ e{p; s') = — 6ss'- The decomposition of A'^ on the 
vectors e^{p; s) reads 

^^W =j:e''{p;s)<p:{z), (2.22) 

s=l 

where (p^i^) are new unknown functions. In terms of these functions, the constraint ( |2.15| ) 
reduces to {(f)^ = (pl T") 

El^'A'] = 0. (2.23) 

s=l 

This is satisfied, in particular, by field configurations of the form (p'Kz) = 0^hs{z) (no 
summation over s), with constant O'l and arbitrary functions hs{z). Indeed, for such 
fields, the color vectors {<t>a{z)} and {4>a{z)} (s = 1, 2, 3) are parallel in color space for 
any s. In this letter we restrict ourselves to such configurations for the color group SU{2) 
(^fabc _ gibc-j^ ^^^ assume, for simplicity, that O^ = 6^, i.e., 

A^iz) = J2^'{p-^a)ha{z)T\ (2.24) 

a=l 

The functions ha{z) (a = 1, 2, 3) satisfy 

as follows from eqs. ( p.l6| ). By contracting these equations with e^{p]a) and using 
eqs. ( ^.171 ), we get 

{u"^ -p^yiii + nlhi + g'^{hl + hl)hi = 0, 
{uj^-p^)h2 + nlh2 + g^hl + hl)h2 = 0, 
{uj^-p^)h3 + nlh^ + g^hl + hl)h3 = 0. (2.26) 

This system admits the following integral of the motion, 

n = lico'-p')j:hl+^- nlihj + hi) + ^-nlhj + ^ {hj hj + hlhl + hi ht) , (2.27) 



which acts as an effective Hamiltonian for a point particle with coordinates ha (the cor- 
responding conjugate momenta being hg). The potential in ( p. 271) is a positive, strictly 



increasing function of the coordinates ha- Accordingly, the energy conservation prevents 
any trajectory {ha{z)} from getting too far away from the origin. 

The limiting case p = 0, cj 7^ 0, of the system (p.26|) , corresponding to global color 
excitations of the plasma, has been investigated in Refs. P, |^. In this limit, VLl = Qt = 
ujp, and the system becomes symmetric with respect to permutations of the functions 
ha- It is then straigthforward to analytically construct particular periodic solutions which 
take advantage of this symmetry. For example, the p = system admits in-phase periodic 
oscillations of the type hi = h2 = h^. More generally, numerical studies of the symmetric 



system [|T0| showed that the oscillations are quasi-periodic for small amplitude, but they 
become unstable as the amplitude is increased. We expect these general properties to 
remain valid for the asymmetric {Qt 7^ ^l) system (p.26|) as well, and, in particular, the 
small amplitude abelian-like plasma waves to be stable with respect to nonlinear effects. 
Note that, since the frequencies Qt and Ql are generally incomensurable, we do not 
expect to find any periodic solution superposing both longitudinal and transverse waves. 
However, periodic solutions of the system (|2.26|) corresponding either to longitudinally, 
or to transversally polarized plane waves can be constructed. 

3 Longitudinal plane waves 

We set hi = h2 = in eqs. ( |2.26|) . We get then a simple harmonic oscillator equation for 

ico^-p'')h + ^lh3 = 0, (3.1) 



which has the general solution h^{z) = Ci cosulz + C2 sinz/^z, with z/^ = ^i/y/uj'^ — p"^ 
and Ci, C2 are integration constants. It is convenient here to replace p^ = {uj,p) by 
k^ = ulp^ = {hP , k), and to set h^iji ■ x) = h^^z). Then, Vl\ = k^k^ and hi satisfies 

hL + hL = 0. (3.2) 



(The dots refer here to derivation with respect to k ■ x.) With these new notations, the 
longitudinal plane- waves that we consider are of the form 

A1{x) = e^(A;; 3) /^^(A; ■ x) T^, (3.3) 

where the frequency fc*^ is related to the wave vector k = |k| by 

k' = 3ul(Q(^) -iV (3.4) 



"K^Kk, 

This last equation, which is easily deduced from ( |2.18| ), is identical to the dispersion 
equation of the abelian-like longitudinal modes|I|]. This is not surprising since the Ansatz 
( |3.3|) led us effectively to a linear, abelian-like, problem. In fact, this property holds for 
any purely longitudinal plane-wave. Indeed, in that case, A'^{z) = e'^{p; 3) (pai^): so that 
[^^,^'^] = for any pair (/i, z/) and the general equations (|2.16|) reduce to harmonic 



equations for each of the functions (p^i^)- The only manifestation of the non-abelian 
structure lies in the constraint ( p.l5| ), which becomes [(f)^,(f)^] = 0. A particular solution 
is 4>a{z) = qahiiz), with constant QaS. For SU{2), this is the only solution, and it is 
equivalent to the Ansatz (|3.3|) , up to a trivial global color rotation (sending the color 
vector QaT"' onto the T^ axis). For a larger color group, other solutions will exist; e.g., for 
SU{3) we can choose (p^iz) = 4>l^{z) T^ + (f)g{z) T^, where T^ and T^ are the commuting 
generators of the group. 

To complete our analysis of the longitudinal modes, we evaluate the corresponding 
energy density and Poynting vector. We return to the particular configuration ( |3.3| ) and 
denote by k'^ = uJiik) the solution of the dispersion equation ( |3.4| )P| [ujiik) > k for any 



/c, ujii^) = Up). The associated field strengths are El{x) = —kJuf{k) — k'^hiT^ and 
Bi = 0. By using these results, together with eqs. (|1.6| )- (|1.7] ), ( p.l2|) -( p.l3|) , and the 
dispersion relation ( ^^ , we obtain 

jT(.)^i{.fit)-e){hl^{-^!L_^2)nl}, (3,5) 



and 



S,ix) = l^\ul-{ufik)-e)\hl\.. (3.6) 



As expected, these quantities satisfy the conservation law 

doTl\x) +d,Sl{x) = 0, (3.7) 

for hi satisfying ( |3.2| ). This last equation also clarifies the physical interpretation of the 
conserved quantity ?i, eq. ( p.27|) . Indeed, since the fields depend on x^ only through k-x, 
it follows that 9o Tf = Wi(fc)Tf and^i^i = -k-S^ = - (A:/cUi(fc))(9o k- S^. Thus, the 



conservation law (|3.7|) shows that the following quantity 



n\x) - ^ k . S,(a:) = \{u^m ' k') [hi + hi) (3.8) 

is constant along the trajectory. This is precisely the integral of the motion 7^^, as 
obtained from eq. (|2.27|) in which we set hi = h2 = and Qf = ujf{k) — k'^. 



4 Transverse plane waves 

We turn now to the more interesting case of purely transverse plane waves, that is, we 
consider the system (|2.26| ) for h^ = 0. The resulting two equations, 

(cj^ - p^) hi + n^hi + g'^ hi hi = 0, 

{uj^-p^)h2 + nlh2 + g^hlh2 = 0, (4.1) 



are symmetric in hi and /12. Again, it is convenient to define ut = ^t/V^"^ "~ P'^ ^^^ 
k^ = ptV^- Then, f2|. = k^k^j_ and from ( p.l9| ) one easily obtains a dispersion equation 
for /c°, 



3 2 



/.2 /i-Oxn 



h.2 _ U2 

rvQ rv 



<t) 



(4.2) 



which is identical to that of the linear (or abelian) transverse modes[0. We denote by 
ujt{k) the positive solution {ujt{k) > k, ujt{0) = Up). The transverse plane waves that we 
are considering here are, in these notations, of the form 

A°. = 0, At{x) = h^^ik ■ x) ei T^ + hl{k ■ x) ea T^ (4.3) 



where ea = ea(k) (recall eq. ( |2.20| )) and h^{k ■ x) = ha{z), for a = 1, 2. The correspond- 
ing field strengths are Et(x) = —ujt{k) Y.a=i,2 h^^aT"-, and Bt(x) = kih'^eiT'^ — 
h}pe2T^j + g h^ph'^liT^, the overdots denoting here derivatives with respect to k ■ x. 
Remark that the three chromoelectric vectors E^ (a = 1, 2, 3) are mutually orthogonal, 
and the same is true for the three chromomagnetic vectors B^ (of course, E3 = 0). Fur- 
thermore, the electric and the magnetic color vectors are orthogonal, E^ ■ Bq = for any 
color a. This is very much similar to the usual transverse electromagnetic waves. 

The energy density and the Poynting vector corresponding to transverse plane-waves 
are obtained by inserting the expressions for the field potentials and strengths above in 
eqs. (pr6D-(pr7D, (|2.12|) - (|2.13|) , and using the dispersion relation (|4^ ): 



T^x) 



{uKk) + k') [{h},f + {h'r) 



.|(. 



1\2 /. 2\2 



T) 



{h: 



T) 



+ 



^m 



?>ujI 



coKk)-k' 



{h^T) + (hr 



2n2 



(4.4) 



and 



Li n2 



L2\2 



St{x) = kujtik)\{h'Ty + {h'Tr + 



One easily verifies that 



u: 



W/ 



UJ„ 



P Vcu|(fc)-P 



-1 



,1 \2 



2\2 



(h'TY + {K) 



(4.5) 



n\x) 



ujt{k) 



k • Srix) = Ht, 



(4.6) 



where the integral of the motion Ht is obtained by setting /13 = in eq. ( p.27|) . 

In order to look for solutions of the system ( |4.1| ), we define h^{k-x) = {Qt / g) fa{k ■ x) 
a = 1, 2, and get the parameter-free system 

,2- 



fi{x) + 

Mx) + 



l+(/2(x)) ]/l(x) 
l+(/l(^))']/2(^) 



0, 
0, 



(4.7) 



for the dimensionless functions /i and /2 . The integral of the motion Tir reads 



Hn 



O'^ 1 



-f^[ft + f^ + f! + fl + f!fi)^-fo 



Q 



To2 
2 



(4.^ 



10 



where Q'^ = uj^{k) — k'^ and ^ is a diinensionless parameter. The Hamiltonian Tix-, to within 
the factor Q^/g'^, is that of a system of two nonhnearly coupled harmonic oscillators with 
coordinates /i and /2 and total energy 6"^. The system ( |4.7| ) has been already analyzed in 
Refs. [^, 1^, where analytic periodic solutions were constructed, as well as in Refs. |T^, 
where the transition from regular to stochastic motion was investigated numerically. 

The simplest solutions correspond to one-dimensional harmonic oscillations along the 
space-color axes 1 or 2. For instance, if /2 = 0, then /i satisfies /i + /i = 0, with the 
general solution /i = ai cos k ■ x + a2 sin k ■ x, and a\ + a\ = 2fP . This is time periodic, 
with period Tq = 2'K/uJt{k). 

Non-linear, but still one- dimensional, periodic solutions are analytically obtained with 
the Ansatz /i = ±/2 = /. These solutions describe in or out of phase oscillations along 
the space-color directions 1 and 2. The function f{k ■ x) satisfies the nonlinear equation 

/ + / + /' = 0. (4.9) 

A particular solution to ( [4.9| ) is 

f{x) = fe cn((202 ^ ^)i/4 f, . ^. ^y (4,10) 

corresponding to the initial conditions /(O) = fg and /(O) = 0. Here, cn(u; n) is the 
Jacobi elliptic cosine of argument u and modulus k, with 

1 / 1 \V2 

1/2 



and the amplitude fg is related to the conserved quantity 9 hj fg = f v2P+T — 1 
The solution ( ^.10| ) is time periodic, with period 

where K{k) is the complete elliptic integral of modulus n. Since | A^l ^ T, 9 ^ 1, and Tg 
remains of order of Tq = 2'K/uJt{k). 
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5 Conclusions 

We have studied here the generahzed Yang-Mills equations in the hot quark-gluon plasma 
for particular plane-wave fields. For the color group SU{2) and for an appropriate Ansatz 
showing correlations between coordinate and color spaces, the problem reduces to an 
effective mechanical system, with three degrees of freedom and nonlinear couplings. Gen- 
erally, this system describes superpositions of longitudinal and transverse plane waves, 
with components along all the three color directions. Simple periodic solutions were ob- 
tained analytically for the particular case where the longitudinal and the transverse plane 
waves decouple. In contrast to what happens at zero temperature, here the plane-wave 
solutions cannot be obtained by simply performing boost transformations on the corre- 
sponding solutions in the comoving framed]. This is so because the mean field equations 
and their solutions are written in the rest frame of the plasma. 
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